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The Biham, Middleton and Levine (BML) model is extended to describe dynamic route 
choices between the residence and workplace in cities. The traffic dynamic in the city 
with a single workplace is studied from the velocity diagram, arrival time probability 
distribution, destination arrival rate and convergence time. The city with double work- 
places is also investigated to compared with a single workplace within the framework of 
four modes of urban growth. The transitional region is found in the velocity diagrams 
where the system undergoes a continuous transition from a moving phase to a completely 
jamming phase. We perform a finite-size scaling analysis of the critical density from a 
statistical point of view and the order parameter of this jamming transition is estimated. 
It is also found that statistical properties of urban traffic are greatly influenced by the 
urban area, workplace area and urban layout. 

Keywords: Cellular automata; BML; residence- workplace; dynamic route choices; origin- 
destination 
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1. Introduction 

In 1992, O. Biham, A. A. Middlcton and D. Levine proposed a two-dimensional 
CA traffic model (BML) to investigate the phase transition and self-organization 
effect in the traffic flow of two dimensions.— The model is defined at a square lattice 
with N x N sites. Within a given time step each site can be held at most by one 
car with three states: moving upward, moving rightward and standing still. The car 
movements are controlled by synchronous traffic lights that allow right moving only 
in even time steps and up moving in odd time steps. At their turn, every car moves 
to its neighbor site ahead whenever the site remains vacant. A sharp jamming 
transition was found between the free phase at low density and the completely 
jamming phase at high density. Recently, many generalizations and extensions of 
the BML model were proposed to investigate various realistic traffic factors in cities, 
such as the anisotropy effect, independent turning and jam-avoiding turning, signal 
control, roads between successive crossings, etc.— ~— 

Moussa incorporated the origin-destination effect of drivers trips in the BML 
model to study the traffic demand in cities.— To characterize the trip, Moussa 
used three different O-D distance probability distribution (ODDPD): exponential, 
uniform, and power-law. Moussa investigated the effect of the ODDPD on the traffic 
flow and the arrival times of travelers. Moussa subsequently studied the evacuation 
processes of cars in cities using above model.— It was found that the ODDPD 
has a great effect on the evacuation processes and the evacuation time diverges in 
obedience to the power law. In-nami et al. introduced a similar model to investigate 
the phase transitions on a decorated square-lattice where the square-lattice point 
and the decorated site denoted the intersection and road, respectively.— In their 
model, a car has a finite deterministic path between the origin and the destination, 
which is assigned to the car from the beginning. 

Huang simulated the city traffic in rush hours in a single-lane roadway with 
open boundaries using the Nagel-Schreckenberg traffic model.— The influence of the 
traffic light control to the traffic flow and travel time was studied. In the morning 
rush hours, downtown parking lots are taken as the destinations for cars moving 
from suburban. When a car arrives at its parking lot, it sinks immediately. In the 
afternoon rush hours, the traffic is in the opposite direction from the morning rush 
hours and original sinks become the sources for cars moving out of the downtown. 
In the present paper, we extend the scenario in road from one dimension to two 
dimensions and introduce the dynamic route choices into car movements within the 
framework of BML traffic model. We find that statistical properties of urban traffic 
are greatly influenced by the urban area, workplace area and urban layout. 

This paper is organized as follows: In Sec. [U we elaborate our model for city 
traffic with the dynamic route choices between the residence and workplace. In 
Sec. El dynamic phase transitions for different sizes of the city and workplace are 
studied numerically. We also investigate the traffic dynamic in the city with double 
workplaces for comparison with a single workplace from the velocity diagram, arrival 
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time probability distribution, destination arrival rate and convergence time. Finally, 
we draw some conclusions in Sec. [4] 

2. Model 

In the original BML traffic model, cars will keep on moving along the vertical and 
horizontal lanes without any target sites if traffic lights permit them to pass. In this 
study, we extend the BML model to describe dynamic route choices between the 
residence and workplace on the traffic in cities. We plan to simulate the scenario 
that office workers drive toward workplaces from their residences every morning. 
The present model includes four modules: 
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Fig. 1. Illustration of the dynamic route 
choices in a square lattice with the periodic 
boundary. There are four cars moving from 
their origins (see Oi O2, O3, O4) to destina- 
tions (see Di D2, D3, D4), which is similar 
to Ref. 11. In our model cars can only move 
rightward or upward. If its destination is at 
the left of the origin, or on the low side (see 
O3 and D3), the car must pull off the lattice 
and then enter the lattice again from the op- 
posite side (see along the dash or solid line 
linking O3 and D3). 



Fig. 2. Illustration of the square lattice 
(city) of size L = 8 with double workplaces 
of size M = 2. The two workplaces, de- 
noted by Wi and W2 in gray areas, are ar- 
ranged uniformly in the city with equal ver- 
tical and horizontal distance (with the peri- 
odic boundary). 



(i) Regional division: In our model, the city is divided into two regions: the residence 
and the workplace. The workplace is a square of the plane lying in the middle 
of the city (See the gray area in Fig.rrjEI around which is the residence (See the 
white area in Fig. [T]). 

a For the periodic boundary used in our model, the position of the workplace has no effect on 
simulation results. Just for the convenience of demonstration, the workplace is put in the center 
of the city. 
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(ii) Initialization of simulation: Each car is assigned a site randomly without over- 
laps from the residence as its origin, and a site randomly from the workplace 
as its destination. The same destination for different cars in the workplace is 
allowed. The origin-destination information is assigned to every car at the be- 
ginning of simulations and all cars start out simultaneously on each way to their 
destinations. 

(iii) Dynamic route choices: In our model, cars select their shortest pathes with a 
simple method similar to Ref. 15 (or see Fig.l). There are two pathes a car can 
select to go to its destination. The car can move upward firstly and turn to right 
when it reaches the site belonging to the lane of the destination site. Then it 
keeps on going until it reaches the destination (see the solid arrow line between 

02 and D2 in Fig. [T]). On the other side, the car may choose another path (see 
the dash arrow line between O2 and D2 in Fig. [IJ. For the periodic boundary, 
the route may go out of bound and then go in from the opposite side (e.g. from 

03 to D3). Every car can select alternatively from the two types of pathes. In 
particular, if the origin and destination site are in the same column (or lane), 
the car has no choice but to move upward (or rightward) (e.g. from Oi to Di or 
from O4 to D4). Counting these cars out, whose origin and destination sites are 
in the same column or lane, half of cars are assigned randomly to move upward 
and the other half to move rightward on departure to avoid the anisotropy effect. 

(iv) Placement of cars after arriving their destinations: We assume that in the work- 
place, there are enough underground garages to accommodate all cars coming 
from the residence. Every site of workplace can hold an infinite number of cars. 
When a car arrives at a destination site, it disappears immediately without 
spending time on parking. 

3. Simulations and Results 

We carry out simulations in a square lattice (city) of size L the workplace of size M. 
We put a number N of cars randomly all over the residence and leave the workplace 
vacant. The density p referred to in this paper is defined as the initial density of 
cars in the residence, which is denoted as p — N/(L 2 — M 2 ). 

We investigate four typical city sizes L = {64, 128, 256, 512} and two types of 
city layouts, i.e. the city with a single workplace (C-l) and the city with double 
workplaces (C-2). We only consider that the two workplaces are distributed evenly 
in the city, which is shown in Fig. [2] To avoid the anisotropy effect, the number of 
cars moving horizontally at the first time step equals that of cars moving vertically. 

The simulation terminates under three criteria: 

(i) If all cars disappear, i.e. every car has arrived at the destination and no one is 
still on the way, the simulation ends spontaneously. 

(ii) If all cars remaining in the city are jammed with traffic and no one could continue 
to move, we terminate the simulation. 

(iii) If the simulation lasts too long and can't still meet the criterion i or ii, it will 
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be terminated forcibly. We make statistics utilizing data generated from the last 
100 time steps. Limited by our computing power, the upper simulation time is 
set to 100,000 steps. In the following text, one time step is referred as one signal 
phase of moving rightward plus one signal phase of moving upward. 

We mainly investigate four typical scenarios about city growing: 

(i) Scenario-1: the area ratio (r) of the workplace to the city remains at 10%. 

(ii) Scenario-2: the area ratio (r) of the workplace to the city remains at 20%. 

(iii) Scenario-3: the area of workplace remains constant during the city expansion. 
We set size M = 20 for the workplace of C-l and M — 14.14 for every workplace 
of C-2. 

(iv) Scenario-4: the workplace narrows down to include only one site, which is the 
extreme case of Scenario-3. 



3.1. Velocity diagrams 

The velocity of each car can be either 1 or 0. The average velocity (v t ) of cars at 
time t is defined as the ratio of moving cars to all cars leaving in the city at time t. 
We mainly investigate the ensemble average velocity (denoted by (v)) by repeated 
many simulations under the same conditions and averaged over the asymptotic 
configurations. If the simulation terminates by the criterion (JTJ) , i.e. all cars arrive at 
their destinations within the given time, (v) is set to one. If only a part of cars arrive 
at their workplaces and the rest get stuck in traffic jam (criterion jn|), (v) is zeroed. 
If the given time (100,000 time steps) expires and some cars are still running, (v) 
is averaged out by taking (v t ) of the last 100 time steps. 




Fig. 3. The polyline shape of average velocity at the end of simulation is influenced by the 
number of samples in every density. In every density, (a) (v) is averaged out from 10 samples; (b) 
(v) is averaged out from 100 samples. We subscribe (or superscribe) each plus sign with a number 
indicating the times of (v) = 1 (or (v) = 0). Every circle is calculated by taking the mean of the 
times of (v) = 1 and (v) = in every density. The surveyed density increases with the step length 
of 0.02. We find in numerous simulations that the fluctuation of curve still cannot be eliminate 
entirely even after adopting 100 samples (not shown in Fig. |3J. 
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Fig. 4. The average velocity ((v)) versus density (p) in the city of (a) C-l in the Scenario-1; (b) 
C-2 in the Scenario-1; (c) C-l in the Scenario-2; (d) C-2 in the Scenario-2; (e) C-l in the Scenario-3; 
(f) C-2 in the Scenario-3, where C-l and C-2 represent the city with a single workplace and double 
workplaces, respectively. The surveyed density increases with the step length 0.02. Every point of 
(v) is drawn by the mean value of ten implements on the corresponding density. 



In order to weaken the random effects of the simulation, ten implements are 
implemented in every density to calculate the average value of (v). Ten samples 
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may not be enough to obtain (v) accurately. Comparing the velocity diagrams of 
different samples in every density in Fig. [3J we find that with the increase of samples, 
the trend of decreasing monotonically is obvious. We judge that the fluctuation and 
the nonmonotonic trend of the polyline in Fig. |3ja) is due to the randomness of car 
initial positions. However, adopting greater than 10 implements in every density is 
beyond our current computing power, especially for the large city size (e.g. L = 512). 
We have to set 10 samples in every density as a compromise. 

Fig.|4]shows the average velocity ((v)) as a function of the density (p) for the four 
different city sizes and two types of urban layouts. The velocity diagrams display 
that (v) decreases from one to zero gradually with the increase of density and 
vanished beyond the critical density p c . There is a transitional region in the velocity 
diagrams where the system undergoes a continuous transition from a moving phase 
((v) =/= 0) to a completely jamming phase ((v) = 0). Comparing different city sizes, 
we find p c keeps decreasing as L increases coinciding with Ref. 1. 

Comparing the velocity diagrams of C-l and C-2 in the Scenario-1 in Fig.JUa)- 
(b), we find that the critical density of C-2 is about 0.05 higher than that of C-l 
except for L = 64. The result indicates that the jam probability of the urban layout 
of double workplaces is less than that of a single workplace especially for large city. 
The curves of C-2 stand closer by one another than that of C-l, which reveals the 
C-2 is less insensitive to the finite size effect than C-l. 

The velocity diagram in the Scenario-2 is shown in Fig. Htc)-(d), which is very 
similar to the Scenario-1. The p c in the Scenario-2 is greater than Scenario-1 and 
the growth amplitude is larger for C-l than C-2. This indicates the road capacity 
gain from the increment of the workplace area is more considerable for C-l than 
C-2. 

In Fig. BJe)-(f), it is obvious that the traffic is terrible for both layouts in huge 
cities. The p c in the Scenario-3 is extremely lower than Scenario-1, 2 when L = 
128, 256, 512 whether for C-l or C-2. Hence it would be unwise to enlarge the city 
while remain the area of workplace or city center constant. When the city extends 
to L = 128, 256, 512, the p c of C-2 is lower than C-l, which is not discovered in the 
Scenario-1, 2. Why this could happen is still unknown. 

3.1.1. The J^th scenario: the workplace narrows down to include one site 

In this section, we will discuss the Scenario-4, i.e. the workplace narrows down to 
include one site, which is the extreme case of Scenario-3. We assume that this site, 
which acts as a "black hole" , can accommodate infinite number of cars. We mainly 
investigate the velocity diagrams and the typical jam configurations. In Fig. [5j we 
plot the velocity diagrams for the single workplace and double workplaces each of 
which has only one site. 

For the city with a single workplace, all cars can arrive at their destinations 
and never get into completely traffic jam. If the workplace is in one piece including 
more than one site, cars going for different goals may be obstructed each other and 
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(a) (b) 

Fig. 5. The velocity diagrams for the (a) a single workplace with only one site; (b) dou- 
ble workplaces each of which has only one site. In Fig. |5(a)| (v) remains one for city sizes 
L = 64, 128, 256, 512 in any density. 
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Fig. 6. The typical configurations observed at several stages of a certain implement in the city 
of size L = 64 with a single workplace: (a) in original state; (b) after 100 iterations; (c) after 400 
iterations; (d) after 1000 iterations; (e) after 2000 iterations; (f) after 3000 iterations. The black 
point at the top right corner of each subfigure stands for the workplace (See the black points 
(64,64). This is just for the convenience of demonstration. For the periodic boundary used in our 
model, the position of the workplace has no effect on simulation results). On the top of each 
subfigure, i, v, N stand for the current number of iteration, current average velocity, cars leaving 
in the city, respectively. 
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the gridlock will form inside the workplace. But when the workplace includes one 
site, this inner obstruction will not happen and all cars will fall into the "black 
hole" even if the city is full of cars (p = 1). Fig. [6] shows the typical configurations 
observed at several stages of a certain implement in the city of size L = 64. Initially, 
all lattice sites are occupied by cars except the workplace. All cars get ready to go 
to the workplace (that black point). The empty area spreads out over time from 
top right corner to lower left corner. The center of the empty area is almost square, 
which is due to the way of route assignment in our model (See Fig. [T]) . 

For the city with double workplaces each of which has one site, the frozen-jam 
will happen quickly even though the density is very low. In our simulation, we 
randomly assign one of the two sites to half of cars as their destinations. Each 
workplace site can't hold the other half of cars which don't belong to it. Since their 
starting points are randomly distributed, cars leaving for different goals will be 
obstructed each other and go into complete jams soon. Fig. [7]shows the typical jam 
configuration observed in a certain implement in the city of size L = 64. 



i = 89, iV = 1179, p = 0.3 




10 20 30 40 50 60 



Fig. 7. The typical jam configuration observed in a certain implement in the city of size L = 64 
with double workplaces, each of which has only one site. (See the black points (17,17) and (49,49)). 
On the top of the figure, i, N, p stand for the current number of iteration, cars leaving in the city, 
initial density of cars respectively. 



3.1.2. The finite-size effects of our model 

The average velocity shown in Figs. 14151 shows strong finite-size effects, similar to 
what is known for other 2d traffic models .^L£ i 14 1 1 7 As defined in subsection 13. 1[ 
the critical density (p c ) is the density above which typical jamming configurations 
form, and thus ensemble average velocity is nearly 0. The relationship between the 
critical density and lattice size is able to be depicted by the exponential model in 
Ref. 19, that is, 



Pc(L) = f3L 



(1) 
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where (3 is a coefficient, while a is the exponent. We use our simulation results to 
verify this hypothesis and obtain the approximate value of a and (3 with the linear 
regression method, which are shown in Fig. [8] and Table. [TJ In order to conduct the 
goodness-of-fit testing, the sample variance (s 2 ) is adopted, i.e. s 2 = Q/(n — 2), 
where Q is the residual sum of squares (RSS) and n is the sample size. The sample 
variance is an unbiased estimation of fit error variance. 
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(b) 



Fig. 8. Log-log plot of p c (L) with the lattice size L = 64, 128, 256, 512 in Scenario-1, 2, 3. (a) the 
city with a single workplace (C-1); (b) the city with double workplaces (C-2). 



Table 1. Estimation of model parameters and fit error vari- 
ance 

Scenario /3i, ai /3-2, 012 s\ s| 

1 1.25, -0.15 1.02, -0.10 0.000063 0.00043 

2 1.31, -0.16 0.99, -0.09 0.00063 0.00041 

3 13.63, -0.70 26.21, -0.87 0.026 0.017 

4 2.77, -0.65 0.048 



In Table. [U j3\ and a\ are the estimation of the coefficient and exponent of model 
[JJwith the values obtained from C-1. Similarly, /3 2 , 012 are the estimation from C-2. 
s\ and s\ are the sample variance corresponding to C-1 and C-2 respectively. In 
particular, in Scenario-4 of C-1 the complete traffic jam will never happen even if 
the city is full of cars (see subsection 13 . 1 . II or Fig. 5(a) I. So in that case the critical 
density dose not exist, which is denoted by dashes (-) in Table. [TJ 

As shown in Fig. [8] and Table. [1] we obtain a good fit to Eq. (JJJ in the Scenario- 
1,2 whether for C-1 or C-2. But a major error will occur when it is applied to 
Scenario-3, 4. From the simulation data available now, it is unclear whether the 



March 8, 2013 0:32 WSPC/INSTRUCTION FILE mypaper 



A Two- Dimensional CA Traffic Model with Dynamic Route Choices Between Residence and Workplace 11 

exponential model ([l} will work, i.e. the critical density and lattice size will remain 
the linear relationship approximately in the logarithmic coordinates, when the area 
ratio (r) of the workplace to the city in Scenario-1,2 increases higher than 20%. This 
is worth carrying out increasingly comprehensive simulations and further research. 

3.2. Arrival time probability distributions 




200 



400 600 800 1000 
Arrival Time 



200 400 600 800 1000 
Arrival Time 



(a) 



(b) 



Fig. 9. Effect of different city sizes on the arrival time probability distributions at relatively low 
density (p = 0.1) in the city with a single workplace in (a) the Scenario-1; (b) the Scenario-3. The 
curve are averaged out by one hundred samples. Because the range of probability of size L = 512 
is one order of magnitude lower than the size L = 256, its diagram is plotted especially as a inset 
within Fig. |9(bJ| for a clearly presentation. These diagrams are plotted in a histogram form with 
the bin granularity of one time step. 



The effect of city sizes on the arrival time probability distribution (ATPD) at 
relatively low density are shown in Fig. [9] We compare the simulation results in 
the Scenario-1 with Scenario-3 considering four city sizes. The distributions of all 
sizes in Fig. 9(a) are unimodal and symmetrical. The position of every curve peak 
equals approximately its corresponding city size. However, the curves shapes of 
large systems in Fig. |9(b)] are asymmetrical with a flat front part (e.g. L — 256 and 
L = 512) and the maximal arrival times are far greater than that in Fig. 9(a) for 
large systems. 

The arrival time probability distributions in different initial densities in the city 
with a single workplace is shown in Fig. [101 We find that with the increase of density, 
the arrival time of maximal probability increases firstly, and then decreases since 
the density reaches the transition region in the velocity diagram. In the low density 
region of free moving phase, all cars can arrive at their destinations in a relatively 
short time (see the curves corresponding A, B in Fig.[T0|). When the density reaches 
the transition region (see the curves corresponding C, D in Fig. [TO]) . only a part 
of cars can arrive at their destinations and their arrival time steps are extended. 
When the density reaches the region of jam phase, only a small percentage of cars 
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(a) (b) 

Fig. 10. The arrival time probability distribution in different initial densities in the city with 
a single workplace (L = 64, M = 20). (a) Six typical densities p = 0.1,0.4,0.54,0.62,0.8,1.0 
are selected to be surveyed, (b) The ATPD curves corresponding the six densities. The curves 
are averaged out by one hundred samples at every surveyed density. When the traffic gets into a 
complete jam in one realization, the arrival time steps of cars leaving in the city are considered 
to be infinite and not concluded in the statistical data. The diagram (b) is plotted in a histogram 
form with the bin granularity of one time step. 



can arrive at their destinations. Even within dozens of time steps, the rest cars 
self-organize into one large cluster and can not go ahead any more (see the curves 
corresponding E, F in Fig. ITU)). The number of arrived cars equals N ■ p(t)dt, 
i.e. the area under the curve multiplied by the total number of cars, where p is the 
arrival time probability, T max is the maximal arrival time. Considering the discrete- 
valued arrival time, the expression should be transformed to N ■ ~^2i=i* Pi> where 
Pi is the arrival time probability of the i-th time step. 

Fig. ITTI shows the arrival time probability distributions of C-l and C-2 in four 
typical densities. We observe that the peak positions of C-2 move to left a tiny bit 
with respect to C-l. But the maximal arrival times of C-2 are larger than those of 
C-l except for p = 0.8 in Fig. 11(d) Comparing the area under probability curves 
of C-l with C-2, we hnd that the number of arrived cars of C-l are almost the same 



as C-2 in the low density region (see Fig. 11(a) and Fig. 11(b)). The number of 
arrived cars are greater for C-2 than C-l in the intermediate density region (see 
Fig. 11(c) ), but less in the high density region (see Fig. 11(d) ). Fig. 11(d) indicates 
that the traffic will be trapped easily and quickly for C-2 than C-l in the region of 
high density. 

The travel time can't be saved in C-2 comparing with C-l, even longer, which 
contradicts the usual experience in reality. We ascribe the experimental phenomena 
to the model definition which permits only two travel directions, i.e. cars can only 
move upward or rightward. If cars can move in four directions as do in real road 
network, the results may be in accordance with actual traffic. The extension of BML 
model including this character is worth further research. 
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Fig. 11. The arrival time probability distribution at the city with a single workplace (L = 64, 
M = 20) is compared with double workplaces (L = 64, M = 14.14). Four typical densities are 
selected to be surveyed corresponding (a) p = 0.3; (b) p = 0.54; (c) p = 0.62; (d) p = 0.8. The 
curves are averaged out by one hundred samples at every surveyed density. When the traffic gets 
into a complete jam in one realization, the arrival time steps of cars left over in the city are 
considered to be infinite and not included in the statistical data. These diagrams are plotted in a 
histogram form with the bin granularity of one time step. 



3.3. Destination arrival rate 

The destination arrival rate is define as the ratio of cars arriving destinations before 
the traffic gets into a complete jam to all cars in the city. If all cars have arrived at 
their terminals finally, the arrival rate is set to one. 

Fig. [12] shows the destination arrival rate as a function of the density for the 
four different city sizes and two types of urban layouts. The trend derived from 
these six diagrams is very similar to the velocity diagrams shown in Fig. |4l As the 
city size enlarges, the curve moves to left and becomes smoother. The traffic in the 
Scenario-3 is worse than Scenario-1,2. It is important to note that the arrival rate 
of C-2 falls faster than C-1 in the transition region except for in the Scenario-3. 
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Fig. 12. The destination arrival rate versus density (p) in the city of (a) C-l in the Scenario-1; 
(b) C-l in the Scenario-2; (c) C-l in the Scenario-3; (d) C-2 in the Scenario-1; (e) C-2 in the 
Scenario-2; (f) C-2 in the Scenario-3. The surveyed density increases with the step length of 0.02. 
The destination arrival rate in every density is averaged out by ten implements. 



3.4. Convergence time 

In this subsection, we measure the convergence time of the simulation by count- 
ing the iterations before all cars have arrived at their destinations. We compare the 
iterations of C-l with that of C-2 in three types of scenarios and four city sizes 
in Fig. [13] We observe the convergence speed of C-l is faster than that of C-2 in 
almost all cases except for the high density region in some instances (see Fig. |13(d)| 



13(e) 13(f) 13(h) I, which contradicts the usual experience in reality. We ascribe the 
simulation results to the model definition which permits only two travel directions, 
just as the interpretation of the simulation in Fig. 1111 The number of iteration rises 
faster in the high density region than the low density region. This is because the 
interaction between cars becomes strong and car movements become rare as the 
density increases. So the ratio of time spent waiting at the crossroads rises rapidly 
with the growth of density. 



4. Conclusions 

In summary, we extend the BML model to describe dynamic route choices between 
the residence and workplace in cities. We plan to simulate the scenario that office 
workers drive toward their workplaces from their residences every morning. 

We firstly investigate the velocity diagrams considering the four different city 
sizes (L = 64, 128, 256, 512) and two types of urban layouts, i.e. the city with a single 
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Fig. 13. Comparison of the convergence time of C-l with C-2 in three types of scenarios: (a-d) 
Scenario-1; (e-h) Scenario-2; (j-1) Scenario-3. Ten implements arc carried out to find the median 
value of iteration in every density. If the traffic gets into a complete jam in an implement, the 
iterations of this implement are set to infinite. When its median value is infinite, the point can not 
be displayed in the figure (This can explain the discontinuity point in Fig. |13(j)| . In particular, 
the curve of C-2 is absent in Fig. |13(1)| This is because the traffic is completely blocked even in 
the density of p = 0.1. 



workplace (C-l) and the city with double workplaces (C-2). The average velocity 
is greatly influenced by the size of city and the area of workplace. Meanwhile, the 
urban layout has some effect on the average velocity. Wc find it would be unwise to 
enlarge the city while remain the workplace area constant. 

We perform a finite-size scaling analysis of the critical density from a statistical 
point of view and the order parameter of this jamming transition is estimated. We 
adopt the exponential model in Ref. 19 to depict the relationship of the critical 
density and lattice size. The parameters of model are estimated with the method of 
linear regression. We obtain a good fit to the model in the Scenario-1, 2 whether for 
C-l or C-2. But a major error will occur when it is applied to Scenario-3, 4. From 
the simulation data available now, it is unclear whether the exponential model ([1]) 
will work when the area ratio (r) in Scenario-1, 2 increases higher than 20%. This is 
worth carrying out increasingly comprehensive simulations and further research. 

The effects of urban size, car density and urban layout on the arrival time 
probability distributions are also studied. The travel time of drivers can't be saved 
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in C-2 compared with C-l, even longer, which contradicts the usual experience in 
reality. We ascribe the simulation results to the model restriction which permits 
only two travel directions. The extension ol BML model by revising this restriction 
is worth further research. 

In addition, we study the destination arrival rate and the convergence time. 
The arrival rate in C-2 falls faster than that in C-l in the transition region except 
in the Scenario-3. We observe the convergence speed in C-l is faster than that in 
C-2 in almost all cases except for the high density region in some instances. The 
convergence time rises faster in the high density region than the low density region. 

In our simulation, we suppose all cars start out simultaneously on each way to 
their destinations. We know that this is not the case in actual traffic. It is possible the 
distribution of car departure time has a great effect on traffic. The departure time 
differences among all cars perhaps can't be ignored, which needs to be investigated 
further. 

In our model, we assume that every site within the workplace can hold an infinite 
number of cars, which is impossible in actual traffic. The next step we will limit the 
capacity of site to accommodate cars and adopt more flexible route assignments, 
which will be closer to actual traffic. 
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